Vortex core states in a minimal two-band model for iron-based superconductors 
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The pairing symmetry is one of the major issues in the study of iron-based superconductors. We 
adopt a minimal two-band tight-binding model with various channels of pairing interaction, and 
derive a set of two-band Bogoliubov-de Gennes (BdG) equations. The BdG equations are imple- 
mented in real space and then solved self-consistently via exact diagonalization. In the uniform 
case, we find that the d x 2_ y 2-w&ve pairing state is most favorable for a nearest-neighbor pairing 
interaction while the s^^-wave pairing state is most favorable for a next-nearest-neighbor pairing 
interaction. The is consistent with that reported by Seo et al. [Phys. Rev. Lett. 101, 206404 
(2008)]. We then proceed to study the local electronic structure around a magnetic vortex core for 
both d x 2_ y 2-wsve and s a .2 y 2-wave pairing symmetry in the mixed state. It is found from the local 
density of states (LDOS) spectra and its spatial variation that the resonance core states near the 
Fermi energy for the d x 2 _ H 2-wave pairing symmetry are bound while those for the s B 2„2-wave pair- 
ing symmetry can evolve from the localized states into extended ones with varying electron filling 
factor. Furthermore, by including an effective exchange interaction, the emergent antiferromagnetic 
spin-density-wave (SDW) order can suppress the resonance core states, which provides one possi- 
ble avenue to understand the absence of resonance peak as revealed by recent scanning tunneling 
microscopy experiment (STM) by Yin et al. [Phys. Rev. Lett. 102, 097002 (2009)]. 

PACS numbers: 74.25.Qt, 74.25.Jb, 74.50.+r, 74.20.Rp 



I. INTRODUCTION 

Recent discovery of the iron-based superconduc- 

t or „1.2,3.4,5.6.7.8.9 .10.11. 12,13.14,15. 16.17.18.19.20,21, 22 nag g en _ 

eratcd considerable interest in the condensed commu- 
nity^ It is the only class of non-cuprate materials with 
high superconducting transition temperature. The iron- 
based superconductors bears both similarity to and im- 
portant differences from cuprate superconductors. On 
the one hand, both families of superconductors have 
a layered structure and exhibit magnetism in undoped 
phase. On the other hand, although the undoped cuprate 
superconductors are antiferromagnetic Mott insulators, 
almost all iron-based parent compounds are still met- 
als, which have an SDW instability with a small mag- 
netic momenti 24 ' 25 In addition, band structure calcula- 
tion a 26 ' 27 : 28 based on the local density approximation 
to density functional theory have emphasized the multi- 
orbital nature in the new superconductors, in contrast 
to the cuprates, where only the d I .2_ y 2-orbital is most 
important. Now it seems commonly accepted from these 
calculations that the Fermi surface consists mainly of two 
electron sheets and two hole sheets in most of the parent 
compounds. 

To uncover the mechanism of superconductivity in 
these materials, the determination of pairing symme- 
try of the superconducting order parameter is a good 
starting point. Up to now, many types of pairing sym- 
metry have been suggested, ranging from the p-wave 
symmetry^ mixed d-w&ve symmetry ; 30 : 31 and extend 
s-wave symmetry ; 32 : 33 besides the ±s-wave symmetry. 34 



Experimentally, the pairing symmetry can be either re- 
vealed by such phase-sensitive methods as multi-crystal 
junctions ; 35 : 36 ' 37 or inferred by non-phase-sensitive tech- 
niques like angled-resolved photoemission spectroscopy^ 
and tunneling^ and other more traditional thermody- 
namical measurements. The first type of technique has 
not been applied successfully to the iron pnictides yet. 
The latter type of techniques are making strides by pro- 
viding detailed information about the nature of single- 
particle excitations. At this moment, the consensus of 
the pairing symmetry in the new class of superconduc- 
tors remains unsettled in view of the fact that conflict- 
ing data have been reported. For example, results from 
some tunneling ; 40 : 41 photoemission^ and nuclear-spin- 
lattice relaxatio n 43 : 44 : 45 : 46 : 47 measurements indicate a d- 
wave pairing symmetry while those from other tunnel- 
ing^ photoemission ; 49 : 50 : 51 specific heat^ 2 - and penetra- 
tion dept h 53 : 54 : 55 measurements show the evidence for an 
s-wave pairing symmetry. 

To identify the pairing symmetry being conventional 
or unconventional, the electronic structure around lo- 
cal inhomogeneity such as single impurities/defects, sur- 
faces/interfaces, and magnetic vortices can provide useful 
information ; 39 : 56 : 57 which can be measured by such local 
probes as STM and nuclear magnetic resonance. Re- 
cently, the local electronic states near impurities 5 ^ and 
surfaces /interface a 59 : 60 '^^ 2 -^ 3 -' 6 ^^ 5 - have been studied in 
the context of s±-w&ve~ or s x 2 y 2-wave 3 ^ pairing symme- 
try in a two-band model for the iron-based superconduc- 
tors, where the existence of impurity resonance or sur- 
face Andreev bound states is discussed. In a single-band 
model as relevant to high-T c cuprates, earlier studie a 73 : 74 
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have revealed the difference of local electronic states in 
the mixed state between s-wave and d-wavc pairing sym- 
metry. The purpose of this paper is to report a com- 
parison study of local electronic structure in the mixed 
state for both the d I .2_ y 2-wave and s^.^-wave pairing 
symmetry within a minimal two-band model of super- 
conductivity. 

The outline of the paper is as follows. In Sec. [Til we in- 
troduce the multi-band model for the iron-pnictide super- 
conductor and derive a set of multi-band BdG equations. 
In Sec. IIII1 we present numerical results first for the uni- 
form case to identify the favorable pairing symmetry for 
a given channel of pairing interaction, and then for the 
mixed state to investigate the local electronic structure 
for both dj.2_ y 2-wave and s x 2 y 2-wave pairing symmetry. 
The effect of the emergent antifcrromagnetic SDW order 
on the resonance core states is also discussed. Finally, a 
summary is given in Sec. IIVI 



By introducing the canonical transformation: 



II. BOGOLIUBOV-DE GENNES EQUATIONS 
FOR MULTI-BAND SUPERCONDUCTORS 

We start with a multi-band model for the normal state 
band structure and introduce the phenomcnological pair- 
ing interactions within each band. The Fermi surface 
topology as mapped out by the magneto-oscillation mea- 
surements consists of weakly corrugated small-size cylin- 
ders. Therefore, we consider only a two-dimensional (2D) 
system for simplification. The magnetic field applied 
along the direction perpendicular to the 2D plane cre- 
ates quantized vortices in the mixed state. The model 
Hamiltonian can be written as: 



ij,a>{3,<7 



(1) 



ij.a/3 



Here the operators Ci aa (c\ aa ) annihilate (create) an elec- 
tron at the i-th site in the orbital a and of the spin 
projection a. In the presence of the magnetic field, 
when the electron hops from the j-site to the z-site, 
a Peierls phase factor is acquired such that tij,a/3 = 
hj,a/3 ex P[* tf^ Jj A - dr], where iy,^ is the zero-field hop- 
ping integral and we assume it to be real, A is the vector 
potential, and the quantity $o = hc/2e is the supercon- 
ducting magnetic flux quantum. In the Hamiltonian (fTj). 
the quantity n is the chemical potential, and Ay iQli g is 
the superconducting pair potential, which is given by 
\j,ap = S a f3Vij ta (ci a iCj a i). In our model, only the spin 
singlet and intra-orbital pairing is considered. However, 
we note that since our model considers one iron atom 
per unit cell, it is consistent with the picture of two elec- 
tron sheets at the M point and two hole sheets at the 
r point obtained from the LDA calculations, where two 
iron atoms per cell are considered. 
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we arrive at a set of multi-band BdG equations: 

= E n 
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Here (u 



f t n ^n^Transpose 



are the eigenstates corresponding 



cigenencrgies E n . The variables with the symbol 
mean matrices or vectors in the orbital space, with the 
single particle Hamiltonian 



H, 



= -t 



ij,af3 



(4) 



and the pair potential subject to the self-consistency con- 
dition 



2 

x tanh 



n(E n >0) 
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(5) 



Notice that the qusiparticle excitation energies are mea- 
sured with respect to the Fermi energy. 

We solve the above set of BdG equations self- 
consistcntly: First guess an initial pair potential Ay and 
exactly diagonalize the equation; use the obtained eigen- 
functions and eigenvalues to update the pair potential; 
repeat the procedure until the desired convergence cri- 
terion is satisfied. Once the self-consistency is achieved, 
the LDOS is calculated as: 



•;„(K)| 2 J(E-E„(K)) 



K,a,n(_E„>0) 

+ \vl a (K)\H(E + E n (K)) 



(6) 



where the factor 2 accounts for the spin degeneracy and 
iV c is the number of magnetic unit cells. This quantity 
is proportional to the differential tunneling conductance 
as measured by STM^ 

The established set of the BdG equations is general 
for describing quasiparticle excitations of a supercon- 
ductor with any number of bands. It is applicable to 
any version of the multi-band tight-binding model as 
recently proposed for the low-energy <i-electron physics 
for the iron-based superconductors i 29 i 67 ' 68 i 69 ' 70 i 71 ' 72 Sev- 
eral group o 29 i 71 ' 72 have pointed out that one may need 
at least three orbitals to accurately reproduce the LDA 
band structure and Fermi surface topology. However, it 
has also been argued recentl y 67 ' 70 that the Fe 3d xz and 
Suborbitals play an important role in the low energy 
physics of these materials. For simplicity, we consider 
in the present pap er a minimal two-band model as sug- 
gested in Ref. [63, with the following tight-binding hop- 
ping integrals t\ = —1.0, t 2 = 1.3, t 3 = f 4 = —0.85. 
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FIG. 1: (Color online) Fermi surface in the Brillouin zone 
with one Fe per cell for various values of chemical potential 
fi = 1.6 (a), 1.8 (b), and 2.0 (c). The diagonal solid lines 
are the nodal lines for d x 2_ y 2-w&ve pairing symmetry while 
the horizontal and vertical dashed lines are the nodal lines for 
s x 2 y 2-wa,ve pairing symmetry. The sheets at the zone center 
and corners are hole pockets, and those at the zone horizontal 
and vertical boundaries are electron pockets. 

These hopping parameters appear in the normal-state 
energy dispersion in the unfolded Brillouin zone as 

p en +£22 //eii-e 2 2\ 2 

Sk = 2 ± V V 2 J +£12£21 ' (7) 

where 

en(k) = — 2ti cos k x — 2i 2 cos k y — cos k x cos k y — fi , 
e 22 (k) = —2^2 cos k x — 2ti cos k y — 4t^ cos k x cos k y — fi , 
ei2 (k) = C2i (k) = — 4*4 sin k x sin k y . 

The chemical potential /i = 1.54 corresponds to the half- 
filled case. In the rigid band approximation, we will re- 
lax this parameter /i to model the superconducting state 
upon the electron doping. The Fermi surface for vari- 
ous values of chemical potential is shown in Fig. [T] As 
expected, the increase of chemical potential, which is 
equivalent to the increase of the electron filling factor, 
enlarges the electron pockets but shrinks the hole pock- 
ets. Throughout the paper, the energies are measured 
in units of |ti|. The temperature is set to be T = 0.01 
for the self-consistency calculations. We further choose 
the identical pairing interaction for each band and denote 
V nn as the nearest-neighbor pairing interaction strength 
and V nnn as the next-nearest-neighbor pairing interaction 
strength. For the calculation of the density of states, the 
Dirac delta-function appearing in Eq. (J5]) is approximated 
by 

with r = 0.01 being chosen. 

III. NUMERICAL RESULTS AND 
DISCUSSIONS 

We first sort out the most favorable pairing symme- 
try in the absence of magnetic field. For the nearest- 
neighbor pairing interaction V nn only, the d x 2_ y 2-wave 
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FIG. 2: (Color online) The bulk density of states as a function 
of energy in a uniform form system with a d x 2_ y 2-w&ve pairing 
symmetry (a) and an s x 2 y 2-wa,ve pairing symmetry for various 
values of chemical potential /i. The pairing interaction for the 
corresponding channel is chosen to be V nn = 1.5 or V nnn = 
1.5. 



(oc cosk x — cos ky) pairing symmetry is more favorable 
than the s x 2 +y 2-wave (oc cosk x +cosk y ) pairing symme- 
try. In detail, no matter whether the initial input order 
parameter is of c? a .2_ J) 2-wave type or s x 2 +y 2-wa,ve type, 
the converged solution will be dominantly of d x 2_ y 2-wave 
type. Similarly, for the next-nearest-neighbor pairing 
interaction V nnn only, the s x 2 y 2-wave (ex cos k x cos k y ) 
pairing symmetry is more favorable over the d xy -wave 
((x sink x sinfc y ) type. This observation is also consistent 
with earlier worki^ 

The typical density of states (DOS) for these types of 
pairing symmetry is shown in Fig. [21 For both types of 
pairing symmetry, the DOS shows a well-defined coherent 
peaks with the position determined by the order param- 
eter. The additional peak outside the coherent peaks is 
due to the van Hove singularity, which is a normal state 
property because its location is fixed in the band struc- 
ture. It shows that the DOS feature in the multi-band 
model can be very different than that in the single band 
model. Specifically, for the d x 2_ y 2--wave pairing symme- 
try, the resultant DOS exhibits a linear behavior within a 
very small energy range near the Fermi energy and then 
reaches a plateau up to a larger energy scale character- 
ized by the amplitude of the c?-wave order parameter. For 
s a .2 1( 2-wave pairing symmetry, the DOS feature is sensi- 
tive to the doping, and can have a y-shape when the 
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FIG. 3: (Color online) The surface and imaging plot of the 
variation of d x 2_ y 2-wave order parameter around the mag- 
netic vortices for band 1 ((a) and (c)) and band 2 ((b) and 
(d)). Here the chemical potential /j, — 1.60, which is corre- 
sponding to the slightly electron doping case. 
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FIG. 5: (Color online) The local density of states as a function 
of energy for the d 2 ,2_j / 2-wave model at at the center of the 
magnetic unit cell (a), at the vortex core center (b), and at 
one of its nearest neighboring sites (c) for various values of 
chemical potential. 



FIG. 4: (Color online) The surface and imaging plot of the 
variation of s x 2 y 2 -wave order parameter around the magnetic 
vortices for band 1 ((a) and (c)) and band 2 ((b) and (d)). 
Here the chemical potential /i = 1.60 is chosen. 



electron pockets are enlarged to cross into the region in 
the Brillouin zone where the s^^-wave order parameter 
has a sign change (see Fig.[5Jb) and Fig.QJc) for fj, = 2.0). 

In the following, we consider the local electronic struc- 
ture around a vortex core for both the d x 2 _ y 2 -wave and 
s a .2 1( 2-wave pairing symmetry in the mixed state. For 
this purpose, the amplitude of the magnetic field is de- 
termined by the condition that each magnetic unit cell 
contains two superconducting flux quanta, i.e., H = 
2&o/(N x N y a 2 ), where N x and N y are the linear dimen- 
sion of a square lattice with lattice constant a. We choose 
N x = 2N y with N y being typically 16 and 20. The num- 
ber of magnetic unit cells is chosen to be N c = M x M y 

N x , or M x = N y /2 and 



N y and M y 



with M x = 
My = N x /2. 

Representative self-consistent results for the order pa- 
rameter of d :E 2_j ) 2-wave and s x 2 y 2-wa,ve pairing symmetry 
are shown in Figs. [3] and [4] The superconducting order 
parameter vanishes at the vortex center and starts to in- 
crease at the scale of superconducting coherence length 
to its bulk value. Since the hopping integrals along the 
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FIG. 6: (Color online) The local density of states as a func- 
tion of energy for the s a .2 !/ 2-wave model at the center of the 
magnetic unit cell (a), at the vortex core center (b), and at 
one of its nearest neighboring sites (c) for various values of 
chemical potential. 
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FIG. 7: (Color online) The imaging of the vortex core states 
at the resonance energy of the peak above the Fermi energy 
shown in Fig. Efb) for the chemical potential values fj. = 1.6 
(a), 1.8 (b) and 2.0 (c). 



FIG. 8: (Color online) The imaging of the vortex core states 
at the resonance energy of the peak above the Fermi energy 
shown in Fig. [SJb) for the chemical potential values fi = 1.6 
(a), 1.8 (b), and 2.0 (c). 



x- and y-dircctions arc noncquivalcnt within each band, 
the order parameter component associated with each in- 
dividual band shows a two-fold symmetry. This is differ- 
ent from the case for the single-band model for high-T c 
cuprates, where the spatial dependence of d-wave order 
parameter shows a four-fold symmetry. However, the 
profile of the order parameter associated with each indi- 
vidual band is related to that of the other band by a 90° 
rotation. 

In Figs. [5] and [51 we show the LDOS at a site mid- 
way between two nearest-neighboring vortices (i.e., the 
center of the magnetic unit cell), at the vortex core 
center, and at one of its nearest neighboring sites for 
both d a .2_ y 2-wave and s x 2 y 2-wavc pairing symmetry The 
LDOS at the midpoint between two nearest-neighboring 
vortices exhibits a grossly similar feature to the case in 



the absence of magnetic field (Compare Figs. [5ja) with 
Fig.^a), and Fig.^a) with Fig.^b)). However, fine os- 
cillating structure is observed in the flat region of density 
of states. By varying the size of the magnetic unit cell, we 
have numerically verified that the oscillation period is in- 
versely proportional to the magnetic field, which suggests 
the fine structure is related to the Landau oscillation. For 
both the d x 2_ y 2-w&ve and s x 2 y 2-wave pairing symmetry, 
two resonant peaks show up near the Fermi energy in the 
LDOS at the vortex core center (see Figs.[S{b) and[5Jb)), 
which reflects the existence of Andreev bound states in 
the limit of isolated vortices. We note that, at the core 
center, the intensity of the peak located at the positive 
energy is much larger than that of the peak located at 
the negative energy When the LDOS is measured near 
the core center (see Figs. [SJc) and^c)), the resonance 
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FIG. 9: (Color online) The spatial variation of the supercon- 
ducting order parameter (a) and the absolute value of the 
SDW order parameter associated with an individual band 
for J nn — 1.95. The other parameter values are fixed at 
Vnnn = 1-5 and fl = 1.6. 



peak below the Fermi energy can be seen more clearly. 
The peaks are shifted toward the Fermi energy as the 
chemical potential is increased. Notably, the peaks in the 
LDOS for the d x 2_ y 2~wave pairing symmetry is not rigor- 
ously located at the Fermi energy, and the peak intensity 
only degrades sightly as the chemical potential varies (see 
Fig. EJb,c)). The peaks in the LDOS for the s x 2 y 2-wave 
pairing symmetry are broadened significantly as they are 
shifted toward the Fermi energy as the chemical potential 
is increased (see Fig. EJb,c)). These results suggest that, 
in the multi-band model, the vortex core states for the 
d x 2_ y 2-wave pairing symmetry arc really bound states, 
while the vortex core states for the s^^-wave pairing 
symmetry can be either localized or extended, depend- 
ing on the chemical potential (i.e., the electron filling 
factor). The evidence for this observation is further en- 
forced by investigating the spatial variation of the LDOS 
at the positive resonance energy, as shown in Figs. [7] and 
[51 where the resonance state for the d x 2_ y 2-w&ve pair- 
ing symmetry is bound around the core center while that 
for the s x 2 y 2-w&ve pairing symmetry begins to show long 
tails as the chemical potential is increased. The results 
are different from the single band modeljS where the 
resonance peak is much broadened in the LDOS for the 
d x 2_ y 2-wa,ve pairing symmetry while it is much sharp for 
the s-wave pairing symmetry, demonstrating the sensi- 
tivity of quasiparticlc properties to the electronic band 
structure. 

Experimentally, recent STM measurements^ on 
BaFei_ a; Co a: As2 observed no zero-bias conductance peak 
at the vortex core center. On the one hand, this exper- 
imental finding is in striking contract to the theoretical 
results above for either d x 2_ y 2-w&ve and s x 2 y 2-w&ve pair- 
ing symmetry. On the other hand, we notice that the 
superconductivity is in close proximity to the magnetism 
in iron pnictides. In particular, the temperature-doping 
electronic phase diagram shows that the antifcrromag- 
netic SDW state is either coexistent with ; 76 ! 77 or incip- 
iently interrupted by^ or disappears immediately be- 



FIG. 10: (Color online) The local density of states as a func- 
tion of energy at the s x 2 y 2-wa.ve vortex core center for various 
values of exchange coupling strength J nn . The other param- 
eter values are fixed at Vnnn = 1.5 and /i = 1.6. 

foreS the superconducting state. In the following, we 
introduce an Heisenbcrg-likc exchange interaction term 
between the electronic spin density to take into account 
the magnetism. We will also restrict the consideration to 
the case of s^^-wave pairing symmetry and assume the 
SDW order comes from the nearest neighbor spin-spin 
interaction. In the mean field approximation, this part 
of Hamiltonian can be written as: 

HsDW = ^ As£)VV,i,a(«iQT — n ial) j (9) 
i . a 

where the SDW order parameter is given by 

&SDW,i,a = 2 J i"+5 m AF,i+S,a > (10) 

s 

with the magnetization niAF.i.a = (nial) ~ ( n iai)- To en- 
force a collincar (71", 0) SDW state, we take the exchange 
coupling along the x- and y-direction to have opposite 
signs, which is consistent with the recent band structure 
calculations on LaOFeAs^ We assume Jx°^ = —Jy = 
Jn°n with J^n to be positive. Although this assumption 
is oversimplified, it should still serve well to demonstrate 
the effect of magnetic ordering. By including this term 
in the BdG equations, we repeat the numerical calcula- 
tions for Varying values Of Jnn (Jnn — Jnn — Jnn)- Li 

Fig. [51 we show a representative spatial variation of su- 
perconducting order parameter and SDW order param- 
eter associated with one of the bands for J nn = 1.95, 
Vnnn = 1.5, and /i = 1.6. The SDW order parameter has 
maximum at the vortex core center and then decreases 
away from the core center, which means the SDW order 
is induced by the magnetic field and is nucleated at the 
vortex core center. In this case, the value of the super- 
conducting order parameter far away from the core center 
is almost unchanged as compared to the case in the ab- 
sence of the magnetism. When the exchange coupling 
is increased, the SDW order becomes more spread while 
the superconducting order parameter becomes more sup- 
pressed accompanied by an expansion of the vortex core 
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in size. We note that for larger values of J nn , the SDW 
can coexist with the superconducting order even in the 
absence of the magnetic field. The magnetic field fur- 
ther suppresses the superconducting order parameter in 
the formation of magnetic vortices. In Fig. 1 101 we show 
the LDOS as a function of energy for various values of 
the exchange coupling strength for the s x 2 y 2 -wave pairing 
symmetry with fixed values of V nnn = 1.5 and fi = 1.6. 
We find that as the SDW order begins to nucleate at the 
vortex core, the resonance peak near the Fermi energy 
in the LDOS is split into a double-peak structure, and 
the overall intensity of the structure is decreased (see the 
curve for J„„ = 1.95). With the increasing exchange 
coupling strength, the SDW order is further enhanced, 
and the peak is further split with decreased intensity 
(see the curve for J nn = 2.00). As the exchange cou- 
pling strength is sufficiently large, the resonance peak 
can be almost completely suppressed (see the curves for 
J nn = 2.05 and 2.10). These results provide one possible 
account for the absence of resonance peak at the vortex 
core center as revealed by the STM on the iron-pnictidc 
superconductor BaFei_ x Coa;As2. 

IV. SUMMARY 

In summary, we have adopted a minimal two-band 
tight-binding model with various channels of pairing in- 
teraction, and derived a set of two-band BdG equations. 
The BdG equations have been implemented in real space 
and then solved self-consistently via exact diagonaliza- 
tion. In the uniform case, it has been found that the 
d x 2_ y 2-w&ve pairing state is most favorable for a nearest- 
neighbor pairing interaction while the s^^-wave pairing 
state is most favorable for a next-nearest-neighbor pair- 



ing interaction. We have also studied the local electronic 
structure around a magnetic vortex core for both d x 2_ y 2- 
wave and s x 2 y 2 -wave pairing symmetry in the mixed 
state. It has been shown from the LDOS spectra and 
its spatial variation that the resonance core states near 
the Fermi energy for the d x 2_ y 2--w&ve pairing symmetry 
are localized while those for the s x 2 y 2-wave pairing sym- 
metry can evolve from the localized states into extended 
ones with varying electron filling factor. Furthermore, 
by including an effective exchange interaction, we have 
shown that the emergent antiferromagnetic SDW order 
can suppress the resonance core states. The emergence 
of the antiferromagnetic SDW states, regardless of being 
field induced or pre-existent provides one possible avenue 
to understand the absence of resonance peak as revealed 
by recent STM experiment. Further studies are still nec- 
essary to finally pin down the origin of the missing reso- 
nance vortex core state in iron-based superconductors. 
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